# Metview Macro

# read the GRIB and BUFR data
t2m  = read("t2m_forecast.grib")
bufr = read("/home/graphics/cgi/metview//course_prep/2015/day_2/processing data/solutions/obs.bufr")

# extract (any parameter) as geopoints
obs = obsfilter(
    output : "geopoints",
    data   : bufr
    )

# compute a new set of geopoints whose locations are from the input
# geopoints, and whose values are from the input field
t2m_at_points = interpolate(t2m, obs)

return t2m_at_points
